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Abstract 

The multiphoton states generated by high-gain spontaneous parametric 
down-conversion (SPDC) in presence of large losses are investigated theo- 
retically and experimentally. The explicit form for the two-photon output 
state has been found to exhibit a Werner structure very resilient to losses 
for any value of the gain parameter, g. The theoretical results are found 
in agreement with the experimental data. The last ones are obtained by 
quantum tomography of the state generated by a high-gain SPDC. 
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I. INTRODUCTION 



Entanglement, the non classical correlation between separated quantum systems, repre- 
sents the physical resource lying at the very basis of quantum information (QI), quantum 
computation and quantum communication. The current proliferation of relevant applica- 
tions of quantum entanglement, ranging from the one-way quantum computation [1] to the 
emerging fields of quantum metrology, lithography, etc. [2] strengthens the need of new 
flexible and reliable techniques to generate entangled states with increasing dimension. 

Entangled photonic qubit pairs generated under suitable phase-matching conditions by 
the spontaneous parametric down conversion (SPDC) process in a nonlinear (NL) crystal 
have been central to many applications ranging from teleportation [3] to various quantum 
key distribution protocols [4]. Today the realization of reliable SPDC sources with large 
efficiency, i.e. with large "brilliance", able to generate entangled pure states is of key interest 
as they largely determine the range of applications of the sophisticated optical methods 
required by modern QI. Recently, high brilliance sources have been realized by using bulk 
[5-7] and periodically poled [8] NL crystals, and fiber coupled SPDC radiation [9]. The 
extension of the attained results to higher dimensional QI qubit carriers could lead to new 
information processing tasks [10] and requires the development of appropriate technological 
and theoretical tools. Significant results in this field are the SPDC generation of four- 
photon entangled states with or without " quantum-injection" [11-16]. At the same time, the 
sophisticated level attained by QI processing of photonic qubits has directed a consistent part 
of the experimental endeavour towards the solution of more technical problems. For instance, 
since the losses due to the photon transmission represent a restriction to the realization of 
several QI protocols, classes of entangled states robust against losses have been investigated 
and a few experimental schemes to generate such states have been proposed [17]. In the 
case of three qubits, entangled states can be classified as "GHZ-type" or as "W-type", the 
latter exhibiting a larger robustness in bipartite entanglement [18]. The W-states can be 
generalized to an arbitrary number of qubits and can be generated efficiently by spontaneous 
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or stimulated emission in a high gain (HG) regime, i.e. when the parametric "gain" is large: 
g >> 0. The first experimental realization of 3 photons W-states by spontaneous emission 
have been recently reported [19,20]. 

In this framework the investigation of the multiphoton state generated in HG condition 
is of fundamental importance, both on a conceptual and practical sides, e.g. for non-locality 
test [21,22] or for QI applications. Recently multiphoton entangled states created by HG 
SPDC have been observed [23]. The aim of the present paper is to investigate relevant 
aspects of the output wavefunction of the HG scheme. In Section II, the density matrix 
representing the 2-photon reduced state arising from HG SPDC after propagation over lossy 
channels is analyzed theoretically. We find that, for any value of g, the resulting 2-photon 
state is a "Werner state" [24], i.e., a superposition of a maximally entangled (singlet) state 
with a fully mixed state. This condition is investigated thoroughly by exploiting the exten- 
sive knowledge available on the Werner states and modern techniques like the "entanglement 
witness" and the various forms of state entropy [25,26]. Resilience of entanglement is the- 
oretically demonstrated for any value of g in the high loss {HL) approximation [27]. In 
Section III, the results of the theory are compared with the corresponding experimental 
data obtained by conventional quantum state tomography (QST). 

II. SPDC IN HIGH LOSSES REGIME: GENERATION OF WERNER STATES 

In this Section, we theoretically analyze the effect of losses on HG SPDC multiphoton 
state. Recently Durkin et al. [27] demonstrated the persistence of some kind of symmetries 
implying entanglement in multiphoton SPDC states in presence of polarization independent 
photon losses. Here we explicitly derive the expression of the SPDC density matrix in 
regime of induced high photon losses through coincidence measurements and demonstrate 
that it corresponds to a Werner state for any value of g. In the present model effects of 
both losses and imperfect detections on the output states are simulated, as usual, by the 
insertion of beam splitters on the two propagation modes k; (i = 1,2) (Fig.l). The results 
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of this study demonstrate the resilience of bipartite entanglement for any value of g. This 
implies that, even in absence of these induced losses, the initial SPDC state is entangled 
for any g, because of the basic impossibility of creating or enhancing the entanglement 
by means of local operations acting on non entangled state [28]. The main motivations 
for the present investigation resides in the experimental entanglement assessment on multi- 
photon HG fields by introducing light absorbing filters on the correlated photons paths. The 
approximate expression of the density matrix also provides an intuitive explanation of the 
behavior of SPDC states in the HG and HL regimes. 

The Hamiltonian of the SPDC process in the interaction picture reads [13,29] 

H = in{a\ H a\ v — a\ v a\ H ) + h.c. (1) 

where a\j represents the creation operators associated with the spatial propagation mode 
kj, with polarization j = {H,V}. H and V stand for horizontal and vertical polarization. 
k is a coupling constant which depends on the crystal nonlinearity and is proportional to 
the amplitude of the pump beam. This Hamiltonian generates a unitary transformation 
U = e~ lI r acting on the input vacuum state \0) 1H |0) iy |0) 2if \0) 2V = |0, 0, 0, 0). The output 
state 



xpOUT^ _ jj | Q ) j s eas iiy obtained in virtue of the disentangling theorem [30]: 

^ oc/T ) = ^Ev^TTr 2 r-) (2) 



n=0 



where 



ip'ty is the n-generated pairs term: 
1 n 

i) n _) = — E (-l) m \n - m) 1H \m) lv \m) 2H \n - m) 2V (3) 

and T = tanhg, C = cosh g. The parameter g = Kt int expresses the NL gain of the 
parametric process, being ti nt the interaction time. The average number of photon generated 
per mode is equal to n — sinh 2 g. 

Let us first consider the contribution ip'l ^ ^™ to the overall density matrix. To inves- 
tigate the propagation over a lossy channel, a beam splitter(5S') with transmittivity r\ for 
any polarization and spatial mode is assumed to simulate the effect of channel losses and 



of detector inefficiencies. Furthermore, perfect detectors with tjqe — 1 measure the output 
state [31] (Fig.l). The symmetry of the entangled state after losses is preserved by assuming 
r) to be mode- and polarization- independent. The contribution to the SPDC state is 
expressed in terms of the BS operators {afj_ /Ar } associated with the input modes {kf^}: 



Vw — / ;t ,(^1H-IN^2V-IN ®lV-IN®2H-IN) n 1 , 0, 0, 0) 

1 yn + ln\ 



(4) 



The BS's couple the input modes {kf^} with the transmitted modes {kf UT } and the 
reflected modes {kf UT }. The output state expression is found by substituting the operators 
|ajj_/iv} with their expressions in term of the operators |S^_ oc/T |, associated with the 
output transmitted modes {kf UT }, and the operators jfe^_ 0[/T j, associated with the output 
reflected modes {kf UT }, through the BS input-output matrix [31]: 

\ 

(5) 



a ij-OUT\ l ) 



\ ( ^ 

a: 



ij-IN\t) 



\ ij 



y/fj iy/l-V 

\Wi-v Vv ) 

v \OUT 

The input state ^" ) evolves into an output state ■?/>-) which is defined over the four 
transmitted modes and the four reflected modes. This state is expressed through the Fock 
states: \tiih, niv, n 2 H, n2v) a ®\ n \H, niv, n 2 H, n 2 y) b where the first term in the tensor product 
represents modes transmitted by the BS and hence detected (a— modes) while the second 



term expresses the reflected modes (b— modes). 

The density matrix a OUT = ip™ ^ (iftL 
easily obtained from the previous expressions: 
out 1 / 1 x ■■ 



of the n-pair term of the SPDC state can be 



E Ml E, ;/ / WA}M(M^}) * 



n + 1 \n\ , 

* \h, h, h, h) a <8>\n - k - h,k - l 2 ,k - l 3 ,n - k - h) b * 
* a (A, A, f3, f4\ ®b (n - h - f 1: h - f 2 ,h - f 3 ,n - h - f 4 \ 



(6) 



with 



A(x, {y k }) 



n \ In — x\ x \ x \ n — x 
K x { yi jUUl V4 , 



^_\Yrf 1+V2+m+V4, {—i\Jl — rj) 2n ~ yi ~ y ' 2 ~ y ' A ~ y4 * (7) 
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n n—x x 

and J2x, yi = J2 J2 J2 ■ Since we are interested in the reduced density matrix p™ 

x=0 2/1,2/4=0 2/3,2/2=0 

defined over the transmitted modes {kf UT }, p™= Tr^ 5 we need to trace 



^oc/r over unc l e tected reflected modes. The result is: 



Pa 



1 



E fcA ft}) I'l. *2, '3, Maa (A, A, A, Al 



n+1 \n! 

*S(n — k — li,n — h — fi)S(k — l 2 , h — A)<H^ ~ h,h — fz)8(n — k — l^n — h — f\ 



(8) 



The final expression for the n-pair contribution to the SPDC density matrix is: 

= E Mi E^.(-l) fe+,l (l " V) 2n S(h, k, l 2 )S(h, k, l 3 )S(h, k, h)S(h, k, U) 

\h, h, h, h) (k — h + l-y, h — k + l 2 , h — k + Z 3 , k — h + U\ 



where S(h, k,p) = C^Q = (V) ( fc ^) and C = ^. 

Up to now we have considered arbitrary, polarization-symmetric losses. In the following 
we make the additional assumption of very large losses, i.e. HL, which greatly simplifies 
our task. Such approximation is expressed by the relation rjn « 1, rfn being the average 
number of photon transmitted by the BS per mode. This condition enables us to take into 
account only the terms of the sum Eq.7 of order < r/ 2 , hence considering only matrix ele- 
ment corresponding to no more than two photons transmitted. As final step, we assume to 
detect one photon on each mode kf UT by a 2-photon coincidence technique. In this way the 
vacuum terms affecting one or both vectors k.f UT are dropped. In summary, this coincidence 
procedure guarantees, by post-selection, that we are dealing only with matrix elements aris- 
ing from the tensor product of the states:{|l, 0, 1, 0) , |1, 0, 0, 1) , |0, 1, 1, 0) , |0, 1, 0, 1)}, which 
correspond to the states {\H) 1 \H) 2 , \H) l \V) 2 , \V) l \H) 2 , \V) l |V) 2 }. The n-pair contribu- 
tion Pp OSt to the SPDC 2-photon density matrix hence reads 



/ 



Plost = >(1 - V) 2n C 



71-1)0 

(l + 2n) -(n + 2) 

-(n + 2) (l + 2n) 

(n-1) 



\ 



(10) 



6 



We note that the above density matrix has the form of a Werner state p w = p + 
^-jj^I, with p = which is a mixture with probability p of the maximally entangled 

state |^_) = 2~ 1/2 (\H) l \V) 2 - \V) 1 \H) 2 ) and of the maximally chaotic state J/4 being I 
the identity operator on the overall Hilbert space. These states are commonly adopted in 
QI, since they model a decoherence process occurring on a singlet state traveling along an 
isotropic noisy channel [32]. 

The complete density matrix for the SPDC output state is obtained by substituting the 

oo 

Pp OSt matrices into the expression p{[ = ^ X (n + l)r 2n ppos*- All the terms Pp OSt sum up 

n=0 

incoherently. Let us explain the latter procedure. In the actual conditions any ^™ 
leads, after the BS action, to two transmitted photons and 2(n — 1) reflected photons. Differ- 
ent n number of input pairs lead to the discard of a different numbers of reflected photons, 
hence any mutual coherence is destroyed after the tracing operation. The normalized density 
matrix turns out to be: 

fl-r, n n n \ 



1-p 

4 














1+p 

4 


_p 

2 








p 

2 


1+p 

4 














1-p 

4 



(11) 



The SPDC density matrix pH, given by the sum of Werner states, is a Werner state 
itself, with singlet weight 

1 



p 



(12) 



2r 2 + i 

with T = (1 — rj) tanhg. In the limit 77 — > 0, Y = tanhg. For large values of g , i.e. for r — > 1, 
and in the hypothesis of very high losses, the singlet weight p > | approaches the minimum 
value |. Since the condition p > \ implies the well known non-separability condition for a 
general Werner state, we have demonstrated for large g the expected high resilience against 
de-coherence of the entangled singlet state [24]. The graph of Fig. 2 shows the behavior of 
singlet weight p as a function of the interaction parameter g. 
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III. EXPERIMENTAL REALIZATION 



The previous theoretical results have been experimentally tested for different values of the 
parametric NL-gain g: Fig.3. The main source was a.Ti : Sa Coherent MIRA mode-locked 
laser further amplified by a Ti : Sa regenerative Coherent Rega 9000 device (A) operating 
with pulse duration 180/s. The amplifier could operate either at a repetition rate 250kHz or 
100kHz leading to an energy per pulse, respectively, of 4yuJ and 8/xJ. The output beam was 
frequency doubled in a UV beam at A p = 397.5nm through a Second Harmonic Generation 
(SHG) process, achieved by focusing the infrared beam into a 1mm thick BBO crystal 
(P — barium — borate), cut for type I phase matching, through a lens with focal length of 
20cm. The nonlinear (NL) crystal was placed at 5cm from the beam waist in order to avoid 
crystal damage and beam spatial distortion. The UV beam then excited a SPDC process 
in a L = 1. 5mm thick BBO NL crystal slab: Fig. 3. The SPDC— generated photons with 
degenerate wavelengths (wl's) A = 2A P = 795nm propagated along the k x and k 2 modes. A 
A/2 waveplate (wp) and a ^ thick BBO were placed on each mode to ensure the accurate 
compensation of all residual birefringence effects coming from the main BBO crystal, cut 
for type II phase matching [5]. In each mode kj, an additional glass plate (G p ) ensured 
a tight balance between the two polarization emission cones of the SPDC process. The 
balancement between the two cones was achieved by suitable tilting of G p in order to vary 
the ratio between the transmittivities for the s— and p— polarized waves. Calibrated neutral 
attenuation filters (At) placed in modes k x and k 2 were adopted to assure the condition of 
high-losses and hence single-photon detection regime. The polarization states analysis was 
carried out through two n— analyzers (T 1 and T 2 in Fig.3) each one consisting of a pair of 
A/4 + A/2 optical waveplates, a polarizing beam splitter (PBS), a single mode fiber coupled 
detector SPCM - AQR1A - FC with an interferential filter of bandwidth A A = 4.5nm 
placed in front of it. The combination of the UV A/2 wp (WPp) and PBSp allowed a fine 
tuning of the UV pump power exciting the NL crystal. 

In a first experiment we estimated the gain value (g) of the optical parametric process 
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and the overall quantum efficiencies of the detection apparatus on both modes. The count 
rates of D x and D 2 and the coincidence rate [D 1: D 2 ] were measured for different values 
of the UV power (Fig. 4). The plots of Figure 4- (a) and 4-(b) clearly show the onset of 
the NL parametric interaction with large g, thus implying the generation of many photon 
pairs. The gain value of the process is obtained by fitting the count rates N { of detector 
Di, dependent of the UV pump power, with the function Ni(g) = R jz^^fyr^ • Here rji 
is the quantum efficiency on mode kj and R is the repetition rate of the pump source 
[23]. The maximal value of gain obtained has been found g max = (1.313 ± 0.002) , which 
leads to a mean photon number per mode n = sinh 2 g max = (2.97 ±0.01). In conclusion 
the maximal total number of generated photon on ki and k 2 modes through the SPDC 
process is M = An = (11.89 ± 0.05). By means of the previous fits we could also estimated 
the overall detection efficiencies (rji) on the k x and k 2 modes which depend on the glass 
attenuation, the fiber coupling, the detection quantum efficiencies: rji = (0.016 ± 0.002) and 
r] 2 = (0.014 ± 0.002). By the previous values we find rjn ~ 0.05. 

The main experimental result of the present work is the full characterization of the 2- 
photon state. We reconstructed the density matrix of the generated 2 — qubit state on k x 
and k 2 modes by adopting the Quantum State Tomography method (QST) [33] . The experi- 
mental density matrix p 1 ^ is obtained by first measuring the 2— photon coincidences [Di, D 2 ] 
for different settings of the QST setup, T x and T 2 , and then by applying a numerical algo- 
rithm to estimate the density matrix. In a low gain condition the SPDC state generated on k x 
and k 2 modes is expected to be in the singlet state l^ - ) = 2 -1 / 2 (\H) ki \V) k2 — \ V) ki \H) k ), 




with excellent agreement between theory and experiment: Figure 5-(d). By increasing g, 
the p elements corresponding to \H) ki \H) k2 ki (H\ k2 (H\ and \V) kl \V) k2 ki (V\ k2 (V\ are no 
longer negligible and the detection of two photon with same polarization is a consequence 
of a multipairs condition: (Figure 5-(c), (b), (a)). The experimental results of the density 
matrices pl^ p for different g— values are in good agreement with the theoretical prediction 
pf[; the mean value of fidelity between the four comparison is T = (0.996 ± 0.002) , where 
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The density matrices p^ p can now be adopted to estimate "singlet weight", "tangle" and 
" linear entropy" of the generated state. The density matrix pw of a Werner state is given by 
the expression (11), as said. The singlet weight (p) can be directly obtained by the matrix 
elements as p = (p^ P ) 22 + (Pci P ) 33 - (Pexp) n - (Pci P ) 44 • Werner states are entangled (p > §) 
or separable (p < |), the extreme conditions being the pure singlet (p = 1) and the totally 
mixed state (p = 0). The tangle is a parameter expressing the degree of entanglement of 
the state, which is defined as r = C 2 , where C is the concurrence of the state [34]; r > 
is a necessary and sufficient condition for a 2 x 2 state to be entangled. Another important 
property for a mixed state is linear entropy (S) , which quantifies the degree of disorder, viz. 
the mixedeness of the system. For a system of dimension 4, it results S = | (1 — Tr (p 2 )). 
In case of a Werner state, we have S w = (1 — p 2 ). For Werner states, "tangle" and "linear 
entropy" are found to be related as follows [35]: 

\ Kl-Sv/T^V) 2 for0<^<§ 
t(S w ) = < (13) 

[ for | < S w < 1 

For each experimental value of g, (S, r) are estimated starting from the experimental density 
matrix. The agreement between experimental results and theoretical predictions are found 
satisfactory: Figure 6-(a). 

An alternative method to establish whether a state is entangled or not is based on 
the concept of entanglement witness. A state p is entangled if and only if there exists a 
Hermitian operator O, a so-called entanglement witness, which has positive expectation value 



Tr 



< 



Opsep > for all separable states p sep and has negative expectation value Tr Op 
on the state p [36-39]. For Werner states p w the method proposed in [25,41] leads to the 
following entanglement-witness operator: 



6 W = \ 



( \H) \H) (H\ (H\ + \V) \V) (V\ (V\ + \D) \D) (D\ (D\ ^ 
v + \F) \F) (F\ (F\ - \L) \R) (L\ (R\ - \R) \L) (R\ (L\ y 



(14) 
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where \D) = (\H) + \ V)) and \F) — ^ (\H) — \ V)) express diagonally polarized single 
photon states, while \L) — A= (\H) + i \ V)) and \R) = A=(\H) —i\V)) express left and right 
circular polarization states. The relationship between the expectation value for a Werner 
state Ww = Tr OwPw and the Werner weight p is found to be 

W w (p) = ^ (15) 



has been estimated 



[7], leading to Ww(p) < for p > |. Experimentally Tr Owpl 
through 8 projective measurements (the 6 projectors appearing in (14) and the operators 
{\H) \V) (H\ (V\ , \V) \H) (V\ (H\} for normalization. [26]). In conclusion, for each g value, 
a point of the Cartesian plane of coordinates (p, W) is obtained: Figure 6-(b). The solid line 
reports the theoretical dependence (15). The comparison demonstrates a good agreement 
between the theoretical prediction and experimental results. 

By the different methods described above the entanglement condition has been found 
to be realized for a value of g up to 1.084 ± 0.002 (Fig.5-(c)), corresponding to an average 
number of photons equal to M = 4n = (6.85 ± 0.03) before losses. For higher values of g 
the presence of bipartite entanglement is degraded by decoherence effects, mostly due to 
imperfect correction of the walk-off effect in the BBO crystal and to time distinguishability 
introduced by the femtosecond pump pulse. 



IV. CONCLUSIONS 

In summary, the present work shows that the multiphoton states generated by SPDC 
exhibit a bipartite entanglement even in the presence of high losses, confirming previous 
analysis [27]. An explicit form has been derived for the output two photon state: a Werner 
state. The theoretical result are found to be in very good agreement with experimental 
data. We believe that the present results could be useful to investigate the resilience of 
entanglement in lossy communication. The present approach can be extended to investigate 
the quantum injected optical parametric amplifier (QIOPA) [42,13,14]. 
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Figure 1. Schematic layout. Inset: losses simulation by beamsplitter. 

Figure 2. Werner weight p versus non-linear parametric gain g. The two-photon density 
matrices p\l are reported for some gain values (g = 0.1, g — l,g — 0.3). 

Figure 3. Experimental setup adopted for multiphoton states generation by means of 
SPDC process and characterization by QST. 
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Figure 4. (a) Count rates [Di] as a function of the UV power (arbitrary unit). The 
continuous line expresses the best fit result, (b) Coincidence rates [D 1 ,D 2 \ as a function of 
the UV power. 

Figure 5. Theoretical p{l (left plot) and experimental (right plot) density matrices 
for different gain values. The experimental density matrices have been reconstructed by 
measuring 16 two qubits observables through the two tomographic setups {T,{\. Each tomo- 
graphic measurement lasted a time t and yielded a maximum twofold counts (cc) for the 
\HV) projection (a) t = 1 sec; cc ~ 9300 (b) t = 2 sec; cc ~ 12000 (c) t = 15 sec; cc ~ 2000 
(d) t = 120 sec; cc ~ 1300. 

Figure 6. (a) The tangle parameter r in function of the Entropy S of the state. Red line; 
theoretical plot (13). (b) Witness parameter W = Tr Owp 1 ^ in function of singlet weight 



p. Red line: theoretical plot (15). 
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